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INCOMPATIBILITY OF TIME-DEPENDENT 
BOGOLIUBOV-DE-GENNES AND GINZBURG-LANDAU 

EQUATIONS 
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AND ROBERT SEIRINGER 

Abstract. We study the time-dependent Bogoliubov-de-Gennes equations for gen¬ 
eric translation-invariant fermionic many-body systems. For initial states that are 
close to thermal equilibrium states at temperatures near the critical temperature, we 
show that the magnitude of the order parameter stays approximately constant in time 
and, in particular, does not follow a time-dependent Ginzburg-Landau equation, 
which is often employed as a phenomenological description and predicts a decay 
of the order parameter in time. The full non-linear structure of the equations is 
necessary to understand this behavior. 


1. Introduction 

The Ginzburg-Landau (GL) model |1] is a paradigm for the phenomenological de¬ 
scription of phase transitions in physical systems. In the static case, its relation to 
the microscopic BGS theory [2] was clarihed by Gor’kov j3] (see also j4] and j5] for 
alternative approaches), and a mathematically rigorous derivation of the GL model 
from BGS theory was recently given in BE]. The present work is concerned with the 
analogous question in the time-dependent case. Arguments for the validity of a time- 
dependent GL equation can be found in the literature, starting with IH1I91EQ] in the 
case of superconductors. These attempts were critically analyzed in m where it was 
argued that the equation can only hold in a gapless regime, however; we refer to [12] 
for a thorough discussion. By applying similar arguments in the study of superfluid 
cold gases, it was stated in |T3l IH] that the non-linear time-dependent GL equation 
applies to such systems for temperatures T slightly above the critical temperature. 

The main message of our present work is that a time-dependent GL equation can¬ 
not be derived from the time-dependent Bogoliubov-de-Gennes (BdG) equations near 
the critical temperature, in contrast to the static case. We shall consider the BdG 
equations for a translation-invariant system, which can be derived in a standard way 
applying the BGS approximation, either to the Heisenberg equations of motion for 
the fermion held operators, or the time-dependent Green’s function [121 EH] • We con¬ 
sider a system which is initially close to a thermal equilibrium state near the critical 
temperature, with non-vanishing order parameter. We then show that, in contrast to 

what would be expected from GL theory, the order parameter does not decay in time. 
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Interestingly, this is a purely non-linear effect. If, instead, we consider the correspond¬ 
ing linear equation, then the solution indeed decays exponentially in time, on a time 
scale that can be calculated via the imaginary part of the corresponding resonance 
pole, which turns out to be proportional to the inverse of the distance to the criti¬ 
cal temperature. The non-linear terms formally contribute a term cubic in the order 
parameter, which indeed resembles a GL type equation; however, small denominators 
close to the Fermi sphere invalidate such formal arguments, and prevent the decay on 
all time scales. 

Our claims are mathematically rigorous and are not based on any ad-hoc approxi¬ 
mations; they are confirmed numerically in [16] in the case of a one-dimensional system 
with contact interactions. 

The present work can be viewed as a continuation of a recent series of studies of 
the mathematical properties of the BCS theory of superconductivity and superfluidity 
[m HE] cni ED]. It is motivated by the current interest concerning the applicability 
of the theory to cold gases, in particular concerning the BCS-BEC crossover [2T]. In 
the BCS regime [221 E3], a rigorous proof of the emergence of a static GL equation 
close to the critical temperature was given in m- In the BEG regime, the prediction 
[H Ea Ea El] of the emergence of the Gross-Pitaevskii equation in the low-density 
limit was rigorously established, both in the static [26] and the dynamical case [27] . 


2. Model and Main Result 


The starting point of our analysis is the BCS model [21 EaE]. The state of a 
(translation-invariant, three-dimensionafl system is described in terms of two quanti¬ 
ties, the momentum distribution 7(k) = (aj^Ok) and the pair density d(k) = (okO^k), 
determining the Cooper pair wave-function via Fourier transform as a(x — y) = 
(27r)“^/^ ^ d^k. They can be conveniently combined to the 2x2 matrix 


F(k) = 



d(k) \ 
1 -7(-k)y 


which satisfies 0 < r(k) ^ for every k. We suppress spin in our notation; the 
pair density a is assumed, for simplicity, to be a spin singlet. The equilibrium state at 
temperature T ^ 0 and chemical potential /i is determined by minimizing the pressure 
functional 

jr(r) = r {e - fi)-f{k)d^k+ ( |a(x)|V(x)d3x-r^(F), (1) 

Jr3 Jr3 


with entropy 

S{T) = - r Trc2 [F(k) InF(k)] d^k. 

Jr3 


^Our analysis extends to one- and two-dimensional systems in a straightforward way. 
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We find it convenient to choose nnits snch that h = 2m = 1, with m denoting the 
particle mass. In ([T]), F(x) denotes a local two-body interaction potential, as appro¬ 
priate for the description of snperflnid gases. Onr results can easily be generalized to 
include non-local interactions as well, which effectively arise via interactions through 
phonons in solids, for instance. 

It was shown in ca that the critical temperature T^. for the model ([1]) is the unique 
T for which the operator KT{—iV) + l^(x) has 0 as its lowest eigenvalue, where 


Krik) 


k'^ — 

tanh j 


That is, for T ^ the pressure functional ([T]) is minimized by the normal state rn(k) 
with 7 n(k) = (l-l-e^*^ and dn(k) = 0, while for T < Tc the pairing density d(k) 

does not vanish identically, showing the transition to a superfluid (or superconducting) 
phase. Strictly speaking, this characterization only applies if Tc is strictly positive, 
which we shall assume henceforth. We shall also assume that + l^(x) has a 

unique normalized eigenvector 0 ^= corresponding to the eigenvalue 0; for radial V^(x) 
this corresponds to the assumption that is an s-wave@ Smoothness and decay 
properties of are analyzed in detail in [6l App. A]. 

For temperatures T slightly below d(k) is proportional to (Tc—to leading 
order in — T. That is, we can write 


d(k) ='0d*(k)-f ^(k) (2) 

with -0 e C of order — T, and ^(k) « V^c — T for small Tc — T. The order 
parameter is, in fact, determined by minimizing the GL type expression 

^gl(i/') = Ggl(T-T,)|i^P + H^ (3) 

for a suitable constant Gql > 0. This follows from the analysis in [6], which is actually 
more general and not restricted to the translation invariant case considered here. 

The time-dependent BCS equation, also known as the BdG equation, has the form 


idtTt{k) = [i/i(k),rt(k)] 


(4) 


with effective Hamiltonian 


Ht{k) 


P-/i At{k)\ 
A,(k) /i-W’ 


where At(k) = 2{27r)~^^‘^ ^V{k — k')at(k')(Tk'. It can be derived from the Heisenberg 
equations of motion for the fermion held operators, applying the same BGS-type ap¬ 
proximations as in the static case [28l [121129] . Alternatively, it also follows from the 
time-dependent Green’s function method introduced in [15] (see also [30]). While only 
certain interaction terms are retained in the BGS approximation, the equation (jl]) al¬ 
lows for pair-creation and annihilation and hence the superhuid density is in general 


^For an analysis of the BCS model without this non-degeneracy assumption, see |20) . 
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not a conserved quantity. Note that Ht depends itself on F* through At, hence this 
equation is non-linear. It is also important to note that the pressure functional ([I]) is 
preserved by the time evolution (0]). 

In this paper we study the time evolution (jll) for initial states close to the normal 
state, for temperatures close to Tc. Closeness is measured in terms of a small param¬ 
eter, which we call h. Let us assume that |T — and also that the initial state 

Fq satishes 

J-(Fo) - J-(Fn) < Ch^ (5) 

for some constant (7 > 0 (independent of /ij. This is satished, for instance, by thermal 
equilibrium states in the temperature range considered, but also by states where the 
order parameter ?/) in ([2]) is modihed by a factor of order one. In addition to the 
assumption Tc > 0, we assume that > 0 and that does not vanish identically on 
the Fermi sphere, i.e., 

sup |d*(k)| > 0, (6) 

which is satished generically. 

With at denoting the pairing density of F* we dehne, for every time t, the complex 
number 

i>t = h~^(a:t:\at). 

For convenience, we multiply by in order for -0 to be of order one. Our main result 
states that, for small h, \'ipt\ remains approximately constant, uniformly in time. 

Theorem 1. Let F* he the solution of with initial state Fq satisfying and 
\T — Tc\ < Then there exists a constant C > 0 such that, for small h, 

Wfjtl^ -^ Ch^^^ (7) 

for all times t. 

We remark that the conditions of the theorem allow to take any value of order 
one. The result states that this value remains constant to leading order in h. This 
holds true even for T > T^, in which case the normal state Fn minimizes the pressure 
functional ([T]). In other words, the order parameter does not tend towards the min¬ 
imum of the GL energy ([3]), as would be expected on the basis of a time-dependent 
GL equation of the form [131IH] 

iddtf} = af) + ( 8 ) 

with a e M, 6 > 0 and d e C with > 0. In fact, for T > one has a > 0, in which 
case the solution to ([8]) goes to zero as t —> oo, in contrast to our main result ([7]). 
Moreover, at T —< 0 (but small), one could for instance start with a state with the 
“wrong” fj, i.e., with a;(k) the equilibrium pairing density multiplied by a complex 

^Throughout the paper, we denote by C generic constants, even if they take different values at 
different places. 



TIME-DEPENDENT BDG EQUATIONS 


5 


number of modulus not equal to one, and our theorem states that this structure will 
be preserved at all times. 

We emphasize that our results do not rule out the validity of the time-dependent 
Ginzburg-Landau equation, in general, which has been successfully employed over 
several decades. What they show, however, is that such an equation cannot be de¬ 
rived from the Bogoliubov-de-Gennes equations, which also appear prominently in the 
physics literature. From the point of view of mathematical physics, it thus remains a 
challenging open problem to unveil the relevant additional physical effects which are 
responsible for the possible emergence of a time-dependent GL equation. 

Proof. The proof of Theorem [1] is divided into three steps. 

Step 1. The hrst step is to show that the energy bound J^(r) — J^(rn) ^ implies 
a decomposition of the form 

7(k) = 7n(k) + T]{k) , d(k) = hipa^fk) + ^(k) 

where = h~^(^a^\a} and where 

||e||2<G/I^ Wvh (9) 

for an appropriate constant G > 0. These bounds follow from the analysis of [6]; we 
shall sketch the main ideas here. From the bound [6l Lemma 1] on the relative entropy 
of F with respect to Fj, we deduce that J^(F) — J^(Fn) can be bounded from below by 

(a\{KT + V)a} + Ji^T(k)(7(k) - ^n{k)fd^k (10) 

2T r 

+ — J Trc 2 [F(k)(l - F(k)) - F,,(k)(l - F„(k))]2 d^k. 

For T ^ Tc, we can bound ^ in the first term. Since is, by definition, the 
non-degenerate zero eigenvector of Kt^ + V, with a spectral gap k > 0, we conclude 
that the first term in flTU]) is bounded from below by + V)f} ^ 

this case. Moreover, the second term is bounded from below by 2T\\j — 7 n ||2 since 
Kt ^ 2T. Hence we obtain the second and third bound in ([9]). Moreover, with 
the aid of the last term in flTOD it is not difficult to show that \'if \ < G, concluding 
the proof of ([9]) for T ^ Tc- The case T < Tc is very similar, using the fact that 
Kt ^ + 2(T — Tc) ^ Kt^ — 2hf instead, and we refer to [6] for the details. 

Step 2. Since the dynamics (jl]) is unitary, the eigenvalues of the 2x2 matrix Ft(k) 
are conserved, for all k e A simple computation shows that they are of the form 
1/2 + Si(k), with 

st(k) = (74(k) - + |dt(k)|2. (11) 

That is, also St(k) is independent of t. 

Step 3. Let Fq be an initial state satisfying ([5]) and let F^ be the solution of (jl]). 
Since the pressure functional is conserved, we also have -F(Ft) — J^(Fn) ^ Ch^ for all t. 
Setting, as above, = /i”^<(Q;*|Dt), di(k) = h'iptC(^{k)+^t(k) and 7 i(k) = 7 n(k)-l-r 7 t(k), 
we hnd that (JID holds for uniformly in t. 
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The equation Si(k)^ = so(k)^ can be written as 

r]t{kf - r]o(kf - {r]t{k) - r]o{k)) tanh 

= |ao(k)P - |at(k)p (12) 

using 1 — 27 n(k) = tanh((A;^ — /i)/2T). We integrate this identity over a thin annulus 
of thickness 6 around the Fermi sphere, denoted by 

= {k e : ||k| — ^/JI\ ^ h}. 

Since tanh((A:^ — ^)/2T) = 0 on the Fermi sphere, 

\ tanh^ < C'^^ 

and hence ([9]), together with the Cauchy-Schwarz inequality, implies that the left side 
of flT^ is bounded, after integration over by C{h‘^ + To estimate the right 

side, we bound 

\\h'il}ta^{k) + 6(k)P - \h'4)oa^{k) + ^o(k)P| 

^ h^|dH=(k)|^ 

- 2h|d*(k)| (|?/)t||^t(k)| + iV'oll^o(k)l) 

-iet(k)r-ieo(k)r. (13) 


The assumption ([6]), together with the continuity of d* (which follows from Prop. 2 
in [6], since the latter implies e L^(M^)), implies that for small <5 



^ C6. 


After integration over the right side of flTdl) is thus bounded from below by 

Ch‘^5 Il'i/'tP — I'^oPl — — Ch^ using again the Cauchy-Schwarz inequality and 

(El). Together with the bound on the left side above, this implies that 

Ch^5\\iJt\^ - ^c{h^ + + ^202^^ _ 


The choice S = h leads to the claim ([7]). 


□ 


3. Comparison with Linear Case 

It is interesting to observe that our main result, namely the fact that |V’t| remains 
approximately constant in time, crucially depends on the nonlinear terms in the time- 
dependent BCS equation. Let us explain this point in more detail. The equation for 
at in (jlD is given by 

idt&tik) = 2(P - p)di(k) -h Ai(k) (1 - 7i(k) - 7i(-k)). 

Writing again 7 i(k) = 7 n(k) -t- 7 t(k) this can be rewritten abstractly as 

idtat = LSat - 2{r]t -h VtWat, 


(14) 
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where S denotes the operator KT^—iV) + h^(x), L is multiplication by 2 — 47n(k) = 
2tanh((/c^ —p)/2T), V is multiplication by V^(x) in x-space and rjt is multiplication by 
i]t(k) in k-space. Consider, for simplicity, the case T > T^, in this case, the operator 
S is positive and the solution to (lT4jl satisfies 

at = S-^/^U{t)S^/^ao + 2i ( - s)S^/\r]s + r]l)Vasds (15) 

Jo 


for t > 0, where 

U{t) = 

is the unitary evolution generated by 

In the second term on the right side of (1T5|1 . we can use (IT^ to express rjs in terms 
of as] this leads to a nonlinear equation for at. Let us neglect for a moment this 
second term, and let us focus on the linear dynamics S~^^'^U{t)S^^'^ao. The spectrum 
of coincides with the one of LS. Its continuous spectrum can easily be seen 

to cover the halfline [— 2p; oo), since LS—2{k‘^ — fi) is relatively compact with respect to 
Moreover, for T = T^, LS has an eigenvalue 0 associated with the eigenvector a^, 
which is embedded in the continuous spectrum. Perturbation theory predicts that for 
T > Tc the zero eigenvalue turns into a complex resonance A, with real and imaginary 
parts of the order T — Tc. 

It is particularly simple to hnd the resonance A of if the potential V is 

ra nk one, i.e., of the form V = —\ip')(ip\ for a. ip e which we assume to be radial 

for simplicity. In this case, Tc is determined by (<^9, = 1 and a* is proportional 

to K^^ip. 

To compute A, we use complex dilation. For 6* e M, we define the unitary operator 
u{9) by 

\u{6)ip^ (x) = 

Alternatively, u{6) = with A = x • k + k • x. Assuming (p to be an analytic vector 
for A, we can extend pe = u{6)p to a strip —f3 < Im0 < 0 below the real axis. In this 
way, we can also define the operators Lg and Sg for all —f3 < Im6* < 0. The resonance 
A then satishes the eigenvalue equation LgSgXe = Aye, which is equivalent to 

Multiplying from the left with (ipg\, we obtain 

2e-2«F-2p + A^" 


1 


Te 
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Note that A vanishes at T = Tc- We can nse implicit differentiation with respect to T 
to expand around T = letting 9^0 afterwards. This yields 


._Tc-T 

2T^ 


p.v. 


|(yc(k)|^ cosh (fk 


l^(k)P 


(P-p)KT,(k) 


d^k — i 




l<^(Vh)P 


to leading order in T — Tc, where the integral in the last factor is understood in the 
principal value (p.v.) sense. 

The fact that Im A < 0 suggests that the corresponding state decays exponentially 
in time. In particular, one would expect from the linear evolution ffT6|) that the order 
parameter satisfies 

l^tl ^ 

to leading order in T — Tc, i.e., it decays on a time scale of the order (T — Tc)“^. Such 
a decay was in fact predicted in [121 Il31 HI]- The meaning of the ss sign here can 
be made precise following the analysis of [21], but the details are not relevant here. 
A comparison with the statement of Theorem [1] shows the importance of the second, 
nonlinear term on the right side of flTHl) for understanding the behavior of a*. Let us 
examine it closer. The function r]t(k) is determined by Eq. (Ej). Away from the Fermi 
sphere, the second term on the left side of flT^ dominates, and we have 


7]t{k) - ?7o(k) 


di(k)p - |do(k)p 
tanh 


(17) 


to leading order, which leads to a cubic equation for the evolution of at. In terms of 
ipt this would even resemble the cubic GL term. On the Fermi sphere the denominator 
on the right side of (16) vanishes, however. As a consequence, rjt is much larger, 
of the order ||di(k)p — |Q;o(k)P|^'^^ according to flT^ . The latter expression equals 
h|d*(k)|||^/'ip — fo leading order. Since d* does not vanish on the Fermi 

sphere, we conclude that \'ipt\ — iV’ol; otherwise, r]t(k) would be too large to satisfy 
(pj). This is exactly the mechanism used in the proof of Theorem [H explaining why 
the nonlinear term in fflSj) plays such an important role in determining the behavior 
of ^i. 
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